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Abstract 



Measurement of the nuclear modification factor (Raa) for heavy-flavor production in heavy-ion collisions tests 
predictions for cold- and hot-nuclear-matter effects. Heavy-flavor production in p+p collisions tests pQCD calcula- 
tions and serves as a reference for understanding heavy-flavor production in heavy-ion collisions. Using the PHENIX 
muon-arm spectrometers, the transverse momentum spectra of inclusive muon candidates are measured for p+p and 
Cu + Cu collisions at -^nn = 200 GeV. After subtracting backgrounds, we obtain the measured invariant yields of 
negative muons from the decay of heavy flavor mesons. For p+p collisions, we measure the charm-production cross 
section integrated over p T and in the rapidity range 1.4 < y < 1.9 to be dcr cS /dy = 0.139 + 0.029 (stat) i^'oss ( s y st ) m b- 
This result is compared to a recent FONLL calculation and to a PHENIX measurement at mid-rapidity. For Cu + Cu 
collisions, we measure the Raa for heavy-flavor muons in three centrality bins for 1 < pj < 4 GeV/c, with suppression 
observed for central collisions. We compare our measurement for central collisions to a recent theoretical prediction. 
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1. Introduction 

Heavy quarks are produced in the early stages of heavy -ion collisions and, therefore, are a probe of the hot dense 
partonic matter produced. Accurate measurements serve as a critical test and constraint for predictions concerning 
the energy loss mechanism of partons moving through such novel nuclear matter. This remains a priority for the field 
of heavy-ion collision physics today. Measurement of heavy quark production in p+p collisions tests perturbative 
Quantum Chromodynamics (pQCD) calculations, in addition to serving as a reference for understanding production 
in heavy ion collisions. 

The nuclear modification factor Raa for Cu + Cu collisions using a p+p reference is defined to be 



where N co n is the average number of nucleon-nucleon collisions in the Cu + Cu collision based on a simple description 
of the nucleus|[T]. In the absence of nuclear effects, Raa is expected to be equal to unity. Values less than unity are 
referred to as "suppression." 

It is important to complement existing measurements for open heavy flavor production at mid-rapidity with mea- 
surements at forward rapidities to fully test theoretical models. Heavy quark production is expected to be modified by 
cold-nuclear matter effects such as shadowing and initial-state energy loss. 



^aa - 
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Figure 1 : Invariant yields of negative muons produced from heavy-flavor meson decay as a function of pj in p+p collisions at V* = 200 GeV 
(open squares) and in Cu + Cu collisions at ^fstm = 200 GeV for three different centrality bins (40-94%, 20-40% and 0-20%), scaled by powers 
of ten for clarity (filled symbols). Solid lines are a fit to the p+p yield, scaled by the corresponding number of binary collisions A'collO 



2. Methodology 

This analysis is based on measurements by the PHENIX Experiment of p+p and Cu + Cu collisions at a center-of- 
mass energy per nucleon-nucleon collision of ysNN = 200 GeV. The south and north PHENIX muon arm spectrom- 
eters cover the pseudorapidity range —1.2 < rj < -2.2 and 1.2 < 77 < 2.4, respectively, over the full azimuth. Each 
muon arm contains a Muon Identifier (MuID) consisting of interleaved planes of steel absorber and Iarocci tubes, with 
the five instrumented "gaps" numbered to 4. 

The double differential invariant yield for muons from the semileptonic decay of heavy flavor mesons is: 

cfiW _ 1 Ni-Nc-N F 
2np T dp T dri 2np T Ap T A7] N evt e^A£ 

where Nj is the total number of muon candidates in the bin that reach the last gap of the MuID (Gap 4) and pass all 
track selection criteria; Nc is the number of tracks corresponding to the irreducible hadronic background estimated 
using a hadron cocktail; Nf is the estimated number of surviving misreconstructed tracks; N evt is the number of 
events; Ae is the detector acceptance and efficiency correction; and e^gJ 1 is the BBC trigger efficiency for events with 
a heavy-flavor muon at forward rapidity. 

A data-driven Monte Carlo "cocktail" is used to predict/simulate the irreducible hadronic background Nc- The 
cocktail consists of a collection of individual simulated primary hadrons which are weighted to match RHIC data. 
Specifically, the transverse momentum spectrum for neutral pions is adjusted to agree with PHENIX central arm 
measurements extrapolated to our rapidity region. Other hadrons are thrown so that they conform to hadron-to-pion 
ratios measured by experiments at RHIC. In order to account for any degradation of reconstructed track quality in 
the Cu + Cu analysis, we embed the hadron cocktail tracks in real Cu + Cu events. The "cocktail" is propagated via 
GEANT through a complete simulation of the PHENIX muon arms. In order to explore sensitivity to the simulated 
propagation through the thick material of the MuID, we produced background estimates using both the FLUKA and 
GEISHA hadron shower codes, with a range of modified hadron-iron interaction cross sections. Tracks which stop 
at MuID Gap 2 or 3 with no hit further downstream (Gap 3 or 4, respectively) have a characteristic distribution of 
longitudinal momentum. Decay muons exhibit a sharp peak due to electromagnetic energy loss, but hadrons exhibit 
a broad shoulder. A highly purified sample of hadrons can be obtained from such tracks by requiring p z > 3 GeV/c. 
The reliance on initial distributions for the hadron cocktail is significantly minimized by subsequent "tuning" which 
further adjusts the hadron weights to match the measurements for hadrons from MuID Gap 2 and Gap 3. 
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Figure 2: cc production cross section as a function of rapidity in p+p collisions at s/s = 200 GeV, measured using semileptonic decay muons at 
forward rapidity (circles) |2| and electrons at mid-rapidity (square) [3]. Curves represent a theoretical prediction (4]|5). 



Simulation studies show that number of decay muons which reach MuID Gap 4 linearly increases with the longi- 
tudinal distance of the vertex from the respective muon arm. Hadrons do not exhibit this dependence. The inclusive 
sample of tracks reaching MuID Gap 4 will exhibit a linear dependence of cIN/j/cIzbbc on £bbc- The weights of the 
hadron cocktail can be adjusted to match this slope for each arm for different transverse momentum bins. 

The number of tracks in the hadron cocktail which satisfy all selection criteria is Nc which we subtract from 
the inclusive yield Ni as indicated in Equation [2] After correcting for acceptance and efficiencies, we obtain the 
measured invariant yield of negatively charged muons from the decay of heavy-flavor mesons in p+p collisions at 
v'snn = 200 GeV and for three centrality bins of Cu + Cu collisions at ysNN = 200 GeV. In order to obtain the 
nuclear modification factor and fully exploit the cancellation of common systematic errors, we use Equation [TJ for 
each bin of transverse momentum and for a common hadron cocktail optimization. This is possible in an experiment 
which measures its own reference data for computation of a nuclear modification factor. 



3. Results 

The measured invariant yields of negatively charged muons from the decay of heavy-flavor mesons in p+p colli- 
sions at = 200 GeV and for three centrality bins of Cu + Cu collisions at y^NN = 200 GeV are shown in Figure W 
scaled by powers of 10 for clarity. The solid lines correspond to a fit to the p+p yield, scaled by the appropriate 
number of binary collisions iV co ii corresponding to the Cu + Cu measurements. In order to estimate the full charm 
cross section for p+p collisions, we express the p+p invariant yield as a differential cross section dcr^-jdy after ex- 
trapolating the measurement to pj — GeV/c using fixed-order-plus-next-to-leading-log (FONLL)||4][5l calculations. 
The charm-production cross section dcr C c/dy is then computed using: 

1 1 dcr M - 

d(r r rldy — • ■ (3) 

BR(c -> p) Cud dy 

where BR(c — > p) is the total charm to muon branching ratio and C//£> is a kinematic correction factor based on 
the FONLL calculation. We obtain a charm-production cross section integrated over p T and in the rapidity range 
1.4 < y < 1.9 of dcr C cldy — 0.139 + 0.029 (stat) *o'ogJ ( svs m b- Figure [2] shows the cc production cross section as a 
function of rapidity using this measurement of muons at forward rapidity [2| compared to a PHENIX measurement 
using electrons at mid-rapidity [3|. Figure [3] shows the transverse momentum distribution of /?aa for negative muons 
from heavy-flavor meson decay in Cu + Cu collisions for three centrality bins [2j. The (red) band in Figure|3]shows a 
recent theoretical prediction |6]|7) fory = 1.65 and pj > 2.5 GeV/c which includes the effects of heavy quark energy 
loss, in-medium heavy meson dissociation, and cold nuclear matter effects such as shadowing and initial state energy 
loss of the incoming partons. 
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Figure 3: Transverse momentum distribution of Saa for negative muons from heavy-flavor meson decay in Cu + Cu collisions at yj'NN = 200 GeV 
for the centrality bins: 40-94% (a), 20-40% (b), and 0-20% (c).(2| The band in (c) is a theoretical prediction from (6][7). 



4. Discussion and Conclusions 

Negative muons from the decay of heavy flavor mesons in p+p collisions at -\fs — 200 GeV are used to measure 
the charm-production cross section integrated over pj and in the rapidity range 1.4 < y < 1.9. The result is dcr C c/dy = 
0.139 ± 0.029 (stat)^gQ5g (syst) mb. This result is compatible with a FONLL calculation within experimental and 
theoretical scale uncertainties. Negative muons from heavy-flavor meson decay have also been measured in Cu + Cu 
collisons at ys NN = 200 GeV for the same rapidity and momentum range, and for 3 bins of centrality. Reasonable 
agreement is observed between the invariant yields for the peripheral (40-90%) and mid-central (20-40%) collisions 
and scaled fits to the invariant yield for p+p collisions; however, a suppression is observed at higher pj for the central 
(0 - 20%) Cu + Cu collisions. This effect is quantified by the reported nuclear modification factors /?aa measured 
in three centrality bins. We observe significant suppression of heavy-flavor muon production in central Cu + Cu 
collisions (0 - 20%). The measurement for central collisions is consistent with a recent theoretical prediction for this 
rapidity and transverse momentum range ||6l |7) . 

The PHENIX Experiment has installed new inner silicon vertex detectors will allow separation of charm and 
bottom contributions in future measurements. Such measurements will constrain predictions and help guide theoretical 
understanding as part of an improving understanding of QCD matter at extreme energy densities. 
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